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Differential rotation occurs in conducting flows in accretion disks and planetary cores. In such 
systems, the magnetorot at ional instability can arise from coupling Lorentz and centrifugal forces to 
cause large radial angular momentum fluxes. We present the first experimental observation of the 
magnetorotational instability. Our system consists of liquid sodium between differentially rotating 
spheres, with an imposed coaxial magnetic field. We characterize the observed patterns, dynamics 
and torque increases, and establish that this instability can occur from a hydrodynamic turbulent 
background. 

PACS numbers: 47.65.+a, 47.20.-k, 95.30.Qd, 91.25.Cw 



Diverse astrophysical and planetary phenomena in- 
volve the close interplay of rotation and magnetic field 
generation. There is theoretical and computational ev- 
idence for a magnetorotational instability 0, Q where 
magnetic fields destabilize differential rotation in elec- 
trically conducting flows 0, IB, 0- This instability is 
conceptually related to centrifugal instabilities such as 
the development of Taylor vortices, but the addition of 
Lorentz forces causes flows to be unstable even when cen- 
trifugally stable. An excellent description of the underly- 
ing process is given in two reviews The instability 
leads to radial outflow of angular momentum, which in 
an astrophysical context implies the enhancement of the 
rate of inward matter flow 0- It also leads to a type 
of magnetohydrodynamic turbulence, different from hy- 
drodynamic turbulence by the forces in balance. The 
nature of and transition to magnetolwdrodynamic tur- 
bulence is of much recent interest The magne- 
torotational instability is thought to affect differentially 
rotating stars and planetary interiors ^|, [lH, and the 
dynamics of accretion disks in protostellar systems and 
around compact objects such as black holes 13]. In a 
geophysical or planetary context, the magnetorotational 
instability could serve as a mechanism for initial field 
growth leading to dynamo states requiring a finite am- 
plitude initiation. 

All our current understanding and appreciation of the 
magnetorotational instability derives from theory and 
numerical simulation P, 0, Until now, there have 

been no direct observations of these phenomena in nature 
or the laboratory, although the suppression of centrifu- 
gal instabilities has been documented 0. In nature, the 
development of this instability occurs in contexts consid- 
erably more complicated than can be simulated, a pri- 
mary problem being the possible presence of small scale 
hydrodynamic turbulence. A central issue we address is 
the growth of the magnetorotational instability within 
turbulent flows, which are unavoidable in experimental 
liquid metal flows in the proper parameter regime. 

The experimental device (see Fig. ^ consists of sodium 



flowing between a rotating inner sphere (radius 
a = 0.050 m), and a stationary outer sphere. The inner 
sphere is made of high conductivity copper, mounted on 
a 0.0125 m radius rotating non-magnetic stainless steel 
shaft which extends along the axis of rotation. The inner 
sphere rotates between 2.5 to 50 revolutions per second 
(= r^o/27r). The 0.010 m thick outer stationary vessel is 
a non-magnetic stainless steel shell (radius 6 = 0.15 m). 
Due to its poor electrical conductivity, the outer vessel is 
relatively passive in the magnetic field dynamics . An 
external magnetic field, from 5^ = to 0.2 T, is applied 
co-axially using a pair of electromagnets. We observe the 
induced fields using an array of 30 external Hall probes 
and the velocity using ultrasound Doppler velocime- 
try 0, '20] . The high electrical conductivity of sodium, 
the highest of any liquid 21], allows significant inter- 
actions between the fluid flow and the currents causing 
induced magnetic fields, and allows an approach toward 
geophysically realistic parameters. 

In addition to the radius ratio a/6, three important di- 
mensionless numbers characterize our experiment. The 
magnetic Prandtl number, Pm = f/rj = 8.8x10"^, is the 
ratio of kinematic viscosity u to magnetic diffusivity rj. A 
small value of Pm is typical of liquid metals, planetary in- 
teriors and other natural systems, but distinguishes them 
from matter in accretion flows, for which Pm ~ 1. The 
Lundquist number, S = Bob/r]{pjii)^^'^ (where p is the 
density and p the magnetic permeability), is the ratio 
of the Alfven frequency to the resistive decay rate. For 
6* ~ 1, system-size magnetic field oscillations within the 
sodium damp in about one period; shorter wavelengths 
are damped more strongly. The magnetic Reynolds num- 
ber, Rm = O^oab/f]^ characterizes the ability of the fluid 
motions to create induced magnetic fields. Our experi- 
ments 22] access the relatively little-explored regime of 
parameter space < 5 < 10 and < Rm < 25. Another 
important (but dependent) parameter is the Reynolds 
number. Re = Rm/Pm- Because of the smallness of 
Pm, the Reynolds numbers for our experiments are large 
(varying between 1.7 x 10^ and 3.4 x 10^), implying well- 
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developed hydrodynamic turbulence. 




FIG. 1: The spherical apparatus and characteristic induced 
magnetic field instability patterns. The device (a) consists of 
a thin stainless spherical outer vessel, a rotating inner copper 
sphere, and liquid sodium in between. The line adjacent to 
V indicates the location of the velocity measurements. The 
resulting dynamics change character under the influence of an 
externally applied magnetic field coaxial with the rotation. As 
the external field increases, one sees states dominated by (b) a 
rotating odd pattern with azimuthal wavenumber m = 1 (01) 
at Rm 26 (50 Hz), 5' = 4.19 (84.6 mT), (c) an even m = 
pattern (EO) at Rm = 10.4 (20 Hz), 5' = 4.37 (88.1 mT), and 
(d) a rotating m = 1 even pattern (El) at Rm = 10.4 (20 
Hz), 5' = 7.08 (14.29 mT). Red and blue indicate outward 
and inward radial (cylindrical) induced magnetic field, and 
green indicates the null values, in an equal area projection of 
the induced magnetic field just outside the outer sphere. 

With no applied magnetic field (our base state), we 
have examined the mean flow along a chord perpendicu- 
lar to the axis of rotation. It has a profile with a velocity 
exponent ( = d\ogQ{r)/d\ogr in the range ~ — 1.4 to 
— 1.6 except in thin boundary layers near the walls (see 
Fig. Here, Q{r) is the rotation rate at a cylindrical 
radius r, and would satisfy Q{r) ~ were ( constant. 
For astrophysical rotation profiles governed by Kepler's 
laws, C = —3/2; that is, Q'^r^ is constant. According to 
the hydrodynamic Rayleigh criterion, flows are linearly 
stable for ( > —2. Nevertheless, we observe 10 — 20% 
turbulent velocity fluctuations in the base state, likely 
generated in the boundary layers, which will be fully de- 
scribed elsewhere. Smaller turbulent fluctuations are also 
observed in the magnetic field in the base state, due to 
interactions between the fluid turbulence and the Earths 
relatively weak field, which is always present in the lab- 
oratory. Profiles with C < are predicted to be magne- 
torotationally unstable, assuming a laminar base state. 
Precise stability boundaries have been calculated for liq- 
uid metals {Pm ~ 10~^) in a cylindrical geometry [13,13 
and for 1 flows in spherical geometry p^. Aside 

from a simple rescaling with P^, the theoretical pre- 
dictions are essentially identical: application of an axial 
magnetic field of strength sufficient to overcome resistiv- 




FIG. 2: The zero field angular momentum and rotation rate 
profiles for Qo/27r = 30 Hz {Rm = 16). From ultrasound ve- 
locimetry measurements, the rising mean angular momentum 
density (a) shows the system is stable to centrifugal instabil- 
ities except for thin boundary layers. The rotation curve (b) 
decreases with cylindrical radial distance from the center, a 
necessary condition for the base state to be unstable to the 
magnetorotational instability. The smooth curve is a Keple- 
rian profile for comparison. The inset (c) shows the velocity 
exponent C with the Keplerian value C = —3/2 indicated by 
a dashed line. 

ity will destabilize long-wavelength magnetorotationally 
driven oscillations. 

Our primary observation is that for fixed rotation rate 
of the inner sphere, above some threshold external mag- 
netic field Bc^ we observe spontaneous excitation of oscil- 
lating magnetic and velocity fields (Fig. These take 
the form of a rotating pattern with azimuthal wavenum- 
ber m = 1 (see Fig. ^Jd). At instability onset the applied 
torque increases (see Fig. 0|, as increased amounts of 
angular momentum are carried from the inner sphere to 
the fixed outer sphere. These observations are evidence 
of the destabilization of the magnetorotational instabil- 
ity. 

The fields and flows consist of toroidal and poloidal 
components |2^. The azimuthal mode numbers m and 
the parity with respect to reflection through the origin 
characterize these components. The primary observed 
magnetic field instability (Fig. ^ is dominated by an 
m = 1 poloidal perturbation that has odd parity. A 
toroidal magnetic field disturbance is also likely, but our 
Hall array, being outside the sodium, can only measure 
poloidal components. We denote the poloidal modes with 
the notation Em (azimuthal wavenumber m with even 
parity) and Om (azimuthal wavenumber m with odd 
parity). After our initial 01 instability, as the external 
magnetic field strength is further increased, a number 
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FIG. 3: Coupled fluctuations in the velocity and magnetic 
fields. Here the amplitude of the induced field for the two 
Gauss modes 01 (top curves) and the amplitude of a velocity 
component (lower black curve) below the outer sphere sub- 
stantiate that Lorentz forces are key to this instability. These 
data were taken at Qo/27r = 7.5 Hz {Rm = 4); the velocity 
was measured 0.115 m from the outer wall. Each mode with 
nonzero azimuthal wave number has two components: when 
they oscillate 7r/2 out of phase the pattern precesses. The 
precession rate depends (nonlinearly) on both applied field 
strength and rotation rate. 



of other modes appear. These include saturated states 
dominated by EO (Fig. [it), El (Fig. [Hi), 02, and E2. 
For applied fields many times Be, the modes often show 
aperiodic changes in pattern, rather than simple preces- 
sion. By varying the rotation rate and external field 
independently, we have navigated the {Rm^S) parame- 
ter plane and determined the regions where these other 
modes dominate (Fig. (S)). 

Our primary instability consists of an m = 1 pattern, 
in contrast to the axisymmetric instabilities that dom- 
inate analogous cylindrical situations. Spherical calcu- 
lations by Kitchatinov and Riidiger have shown m = 1 
instabilities in addition to axisymmetric instabilities jl^ . 
These differences can be understood by contrasting the 
different symmetries of the base states. In the cylindrical 
case, the base state is unchanged by axial translations 
and rotations. Such situations generically show insta- 
bilities to axially periodic patterns, and are known to 
do so for the magnetorotational instability 
Our base state lacks the axial translation symmetry, but 
has approximate rotational and reflectional symmetry. 
In such situations, instabilities involving rotating non- 
axisymmetric patterns (Hopf bifurcations) are generic 
[2^. Indeed, the amplitude of our disturbance shows 
the characteristic rise of a Hopf bifurcation, for both the 
induced magnetic (see Fig. and velocity fields. How- 
ever, the onset is made imperfect by background turbu- 
lence and possibly geometric imperfections. In a system 
similar to ours, Hollerbach and Skinner [2^ found, nu- 
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FIG. 4: The torque increase and variance of induced field 
coefficients as the applied field is varied. We characterize 
the saturated state by the mode with the largest variance. 
These measurements were taken for a fixed rotation rate of 
30 Hz [Rm — 16). The field data reflect fluctuations in the 
coefficients from their means, except the EO state which shows 
significant mean amplitude (squared average shown in green, 
variance of fluctuations in blue), and does so at the same 
applied field values where it shows ample fluctuations. At 
zero field, the base state torque (including confounding errors 
due to seals) here is 1.14 N-m, the black curve shows the 
increase in torque above this base value. 

merically, rotating non- axisymmetric patterns, though at 
much lower Reynolds number. 

Also shown in Fig. [SI are predicted magnetorotational 
stability boundaries, as obtained from an inviscid dis- 
persion relation |23, adapted to the present config- 
uration. The two curves represent the stability bound- 
aries of modes with kb = 1 (red) and kb = 2 (blue) for 
wavenumber /c, i.e. wavelengths of one or one-half a cir- 
cumference. We take the peak rotation rate Qp in the 
dispersion relation to be an estimated fluid rotation just 
outside the inner boundary layer Qp = Qq/^- The cor- 
relation between the predicted stability boundaries and 
the experimental observation of the magnetic field oscilla- 
tions with progressively more spatial structure is strong 
evidence that we have observed the magnetorotational 
instability in our experiment. 

These measurements, beyond being the first direct ob- 
servation of the magnetorotational instability in a phys- 
ical system, have important implications for our under- 
standing of magnetohydrodynamic instabilities. They es- 
tablish that the magnetorotational instability, and a sig- 
nificant increase in angular momentum transport, occur 
in the presence of pre-existing hydrodynamic turbulence. 
We quantify the nonlinear saturated amplitude of the 
angular momentum transport and the patterns and sat- 
urated values of magnetic field over a range of parameter 
values not computationally accessible. Finally, this ge- 
ometry may well be capable, at high Rm, of showing a 
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FIG. 5: Phase diagram of saturated states. Regions are de- 
fined by the mode with the largest variance of fluctuations 
see Fig. m. Some secondary instabilities show hysteresis; 
these data are for increasing S, for fixed Rm- The states have 
regions associated with background turbulence (tb), mode 
01 dominated, EO dominated, followed by 01, El, 02, and 
E2 modes. The lowest Rm and S numbers for these states 
(□) are obtained by extrapolating the maximum saturated 
state amplitude versus Rm (a Unear trend) to zero amplitude. 
Also shown are theoretical stability boundaries for the longest 
wavelength (red) and second longest wavelength (blue) insta- 
bilities, calculated from the magnetorotational dispersion re- 
lation. 

dynamo instability. While the dynamo instability is dis- 
tinct from the magnetorotational instability, it would be 
important to examine a system capable of showing both. 
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